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We give the optimal decomposition of a universal two-qubit circuit using Heisenberg exchange 
interactions and single qubit rotations. Tuning the strength and duration of Heisenberg exchange 
allows one to implement (SWAP) 01 gates. Our optimal circuit is constructed from only three 
(SWAP) 01 gates and six single-qubit gates. We show that three (SWAP) " gates are not only 
sufficient, but necessary. Since six single-qubit gates are already known to be necessary, our imple- 
mentation is optimal in gate count. 
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Introduction In several general solid-state quantum 
computation appoaches0, 0,0, CJ El two-qubit inter- 
actions are generated by a tunable exchange interaction. 
For example, Heisenberg exchange between two electron 
spin qubits results in a (SWAP) " gate, where the expo- 
nent a is controlled by adjusting the strength and du- 
ration of Heisenberg exchange. Single-qubit rotations 
have also been proposed for solid-state computation; no- 
table mechanisms for rotating spin qubits being g-tensor 
resonance 0, @ and localized magnetic resonance 0. In 
general, it is desirable to optimize quantum circuits with 
respect to the number of physical operations required, 
which for most solid-state quantum computation propos- 
als implies that circuits should be optimized with respect 
to the number of (SWAP) 01 operations and single-qubit 
rotations. 

The problem of optimizing quantum circuits for ex- 
ecuting general n-qubit operations is computationally 
intractable. Hence, just as in the classical computa- 
tion case, one needs to develop techniques for optimiz- 
ing only few-qubit circuits and then assemble these cir- 
cuits together in a modular fashion. Toward this end, 
circuit optimization results have mostly dealt with the 
case where controllcd-NOT (CNOT) gates and single- 
qubit rotations are the basic building blocks. For ex- 
ample, for a general two-qubit unitary operation, it has 
been recently shown that three CNOT gates and addi- 
tional single-qubit rotations are sufficient and necessary 
[ToL HH]. Now, it is known that one CNOT gate can be 
realized by two (SWAP) 1 / 2 gates and single-qubit uni- 
tary gatesQ. Hence, six (SWAP) 1 / 2 gates are sufficient 
to implement any two-qubit operation. The question is 
whether this strategy of simple substitution is optimal? 
Or, are (SWAP) 01 gates just as efficient as CNOT gates 
(in terms of gate count) in performing two-qubit oper- 
ations? Our answer to the latter is in the affirmative: 
The (SWAP) 01 gates and CNOT gates are both equally 
efficient at realizing any two-qubit quantum operation 
(when measured in terms of number of gates), and that, 
in order to achieve optimal realizations each type of cir- 
cuit requires its own optimization scheme. 

The primary results of our paper are as follows. An 



arbitrary two-qubit operation can be implemented us- 
ing only three (SWAP) 01 gates and six single-qubit ro- 
tations. We augment this result with a number of lower 
bounds. First, we show that, by considering entangle- 
ment power alone, a CNOT gate requires at least two 
(SWAP) a gates. Next, we prove that three (SWAP)" 
gates are not only sufficient, but in fact necessary, to 
implement an arbitrary two-qubit operation. Our uni- 
versal two-qubit circuit is optimal in the number of both 
(SWAP) 01 and single-qubit gates. 

Heisenberg interaction Let's first fix some notation; 
the four Bell states are written as, l^*) = i(|00)±|ll)), 
|*±) = -^=(|01) ± 1 10>). The SWAP gate is defined as 
SW AP\ip)\4>) = \(f>)\i}). For a C 2 <g> C 2 system, it can be 
written explicitly as, 



SWAP = |$+)($+| + 

+e i7r | *-)(*-|. 
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In this paper, our basic two-qubit gate is (SWAP) a , it 
can be written as, 
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The Hamiltonian of the isotropic Heisenberg exchange 
interaction between electron spins Si and S2 is, 

H = J(t)Si ■ S 2 , (3) 

where S — {o~ Xl o~ y , cr z } is a vector of Pauli matrices, 
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The coupling constant J(t) can in principal be tuned for 
confined electrons[l|. The unitary operator generated by 
this Hamiltonian is, 



exp 



- — Si ■ S2 
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By adjusting the integrated coupling J J(t)dt, the 
unitary operator U12 can naturally realize the gate 
(SWAP) a where a = J J(t)dt/h. In this paper, we 
will use the (SWAP) a gate as the two-qubit gate. It 
was proposed to use the Heisenberg interaction alone to 
implement quantum computing0,|a,@- This scheme en- 
codes one logical qubit as three physical qubits. Addi- 
tionally, one CNOT gate requires 19 Heisenberg inter- 
actions amongst the six physical qubits. We consider 
the scheme where both Heisenberg interaction, as well as, 
single-qubit rotations are available. 

CNOT gate requires two (SW AP) a gates We know 
that one CNOT gate can be realized by two square root 
of SWAP gates, V SWAP. If we use a more general 
gate, (SWAP) 01 , can we realize the CNOT by only one 
(SWAP)" gate and a certain number of single-qubit ro- 
tations? By studying the non-local invariants of the 
quantum gates, Makhlin showed that the CNOT gate 
cannot be constructed by applying the Heisenberg inter- 
action Ji only once, i.e., two (SWAP) a gates are neces- 
sary to construct one CNOT 12] . In this section, we give 
a different proof, using entanglement power, to show that 
the CNOT gate requires at least two (SWAP)" gates. 

The entanglement power of a unitary operator U G 
SU(4:) is defined as, 

E P [U) = Average^ l)8 | V , 2) [S(C/|Vi) ® (6) 

where the average is over all product states \tpi) <8> \^2) G 
C 2 ® C 2 in uniform distribution, see 0], and E is the 
linear entropy which is also the concurrence^). Note 
that for arbitrary Ui ® U 2 G SU(2) ® SU(2), E p (U) = 
E p (Ui U2U). So, the entanglement power of (u\ ® 
vi)(SW AP) a {u2®V2) is actually the entanglement power 
of (SWAP)". 

A simple formula can be used to calculate the entan- 
glement power [ill EH, 

E p (U) = ^-±[(U® 2 ,T l!3 U® 2 T lj3 ) 
+ ({SWAP ■ U)® 2 ,T h3 (SWAP ■ U)® 2 T 1>3 )], (7) 

where Ti j3 acting on C 2 <£> C 2 <S> C 2 ® C 2 is the transposi- 
tion operator: Ti l3 \a, b, c, d) = |c, b, a,c). By tedious but 
straight forward calculations we can show that, 

E p ((SWAP) a ) = 1 - 1 cos(27ra). (8) 

For detailed calculations, see Appendix A in Ref . flif . 
When a — 1/2, (SWAP) a has a maximum entangle- 
ment power of 1/6. Since the entanglement power of 
CNOT is 2/9 [13, which is strictly larger than 1/6, one 
(SWAP) a operator with the help of single-qubit gates is 
not sufficient to realize the CNOT. Hence, at least two 
(SWAP) a gates are necessary to realize a general SU(4) 
operator. 



General two-qubit operation Kraus and Cirac gave 
the following result, (see also an arbitrary unitary 

transformation U G SU (4) has the decomposition, 

U ={u' i ®v' i )e- lH {u 1 ®v 1 ), (9) 

where u\,v\, u' 4 ,v' 4 G SU(2), and 

H = h x a x ®a x + h y a y <g> a y + h z a z ® cr z , (10) 

where ?r/4 > h x > h y > h z > 0. Then g in lflH| can be 
written as, 

H = Aoo|$ + }(* + |+A i|* + )(* + |+A 10 |$-}($-| 

+A u |*-)(*-|, (11) 

with, 

Aoo = hx — hy + h z , A01 = h x + h y — h z , 
A10 = —h x + h y + h z , An = —h x - h y — h z . (12) 

The diagonal form of H thus gives, 

e -iH = e -^oo|$+^$+| + e - iAol |* + >(* + | 

+e- iXl0 \<£-)($-\ + e- iXll \y-)(y-\. (13) 

Vidal and Dawson lQj, Vatan and Williams [llj have 
shown that the operator e~ lH can be realized by only 
three CNOT gates and some single-qubit rotation gates. 
Thus an arbitrary U G SU(A) can be realized by three 
CNOT gates and additional single-qubit gates, see also 
the appendix B in Ref. [18j . 

Arbitrary two-qubit unitary operations require only 
three (SW AP) a gates and six single-qubit gates 

Recall that a CNOT gate can be realized by two 
(SWAP) 1 / 2 gates and a few extra single-qubit opera- 
tions. We know the optimal circuit for a general U G 
SU(4) needs three CNOT gates, so six (SWAP) a gates 
are needed if we simply substitute SWAP circuits for 
CNOT gates. Our aim is to find a circuit to realize 
U G SU(A) optimal in the number of (SWAP) a gates. 

From the result of Kraus and Cirac (Eq. I13[) |15| , we 
need to create arbitrary phases on four Bell states by 
(SWAP) 01 gates. But we already know that (SWAP) a 
applies a phase to the Bell state while leaving the 

other three Bell states invariant. Also, by Pauli rotations 
on one particle of the bipartite state, we can transform 
the Bell states amongst each other. Thus, it is straight- 
forward to create four independent phases on the Bell 
states. We can rewrite the operator exp(— iH) as, 

e -iH = e i{h z -h x -h v )n^+\H^+\ _|_ e 2i(/i„+/i y )|^r-w^-| 
+ e 2i(ft„-h a )|$+\/$+| _|_ e 2l (' l --' l z) ]$-)($- |)(14) 

This operator can be constructed by (SWAP) a operators 
as 

e - iH = e i{hz-h x -h y ) x 
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(7 <g> cT x( r z )(SWAP) 2 ^-^)/ w (/ <g> <r z < 

(I® ^(SWAP) 2 ^*"' 1 ^!® t^) X 

(SWAP) 2{ - hx+hy)l * 



(15) 



This circuit just involves three (SWAP) a operators and 
single-qubit gates which are Pauli matrices. So as a 
whole, we can construct any U £ SU(4) by only three 
SWAP gates and single-qubit rotations. Note that be- 
sides Pauli matrices, general single-qubit rotations are 
also necessary to transform e~ zH to U. 

Up to an overall phase, we can rewrite e~ lH , 

e- lH = (a z ®<j x )(SWAPy(a z ®I) 

(SWAPf(I®<j x )(SWAP) a , (16) 

where a = 2(h x + h y )/ir, (3 = 2(h x — h z )/%, 7 = 2(h y — 
h z )/n. The corresponding circuit for U is illustrated in 
Fig. m 
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FIG. 1: Circuit for arbitrary unitary transformation U £ 
517(4) as decomposed in Eqs. © and (T§). Three SWAP 
gates, and 6 local unitaries (upon combination of u'4 and v' 4 
with Z and X, respectively) are required. 

In any implementation, single-qubit rotations as well 
as two-qubit operations will require physical resources 
such as time and hardware. Hence, it is helpful to con- 
sider the number of single-qubit gates involved in each 
circuit as well. In the circuit of Vidal and Dawson[Io|. 
eight single-qubit gates are used to construct the gen- 
eral U £ SU(4), while in our circuit , six single-qubit 
gates are used. It is known that this is the least possible 
number of single-qubit rot at ions [l^. If we assume that 
each single-qubit rotation is as expensive as a two-qubit 
operation, then our circuit is potentially cheaper. 

We next show that our circuit is optimal in the num- 
ber of {SWAP) a gates used, i.e., three (SWAP) 01 gates 
are necessary to construct a general circuit. Let's write 
out a general unitary operator which contains just two 
{SWAP) a gates, 



U = (U 1 (g>V 1 )(SWAP) a (U 2 ®V 2 ) 
{SWAPf(U 3 ®V 3 ), 



(17) 



where Uj, Vj,j = 1, 2, 3 are single-qubit operations. 

One can group single-qubit unitaries about the 
(SWAP) a operators as follows, 

U = (t/i <g> Vx){SWAP) a (U\ <g> V±) 

(U 2 ® V 2 ){SWAPf{ill ® v2)(U a ® V 3 ), (18) 



where U 2 = IW2, V 2 = V1V2, U 3 = UiU 2 U 3 , and 
V 3 = V\V 2 V 3 . Here we notice that operators (U\ ® 
V x ){SWAP) a {u{®v£) and {U 2 ®V 2 ){SW APf {u\®V^) 
are just SWAP gates in some different basis. So we can 
write this relation as, 



U = (SWAP) a (SWAPf(u <g> v) 



(19) 



A single (SWAP) a gate can create one phase in one 
maximally entangled state, with the orthogonal three di- 
mensional space left invariant. So for two (SWAP) a 
gates, there exist two maximally entangled orthogonal 
states, \xi) and \X2), simultaneously orthogonal to \^~) 
and l^ - ). ^,From the symmetry of the (SWAP) a opera- 
tion, it follows that for every two- (SWAP) "-gate circuit 
there exists at least two orthogonal and maximally entan- 
gled states such that they cannot be assigned a relative 
phase by the circuit. That is, the unitary operator corre- 
sponding to the two- ( SWAP) Q -gate circuit must satisfy 
U\Xj) = ( u ® w )IXj):i = 1>2. Note that local unitary 
operations cannot add a relative phase to two maximally 
entangled states. However, since a general two-qubit op- 
eration can assign independent phases to three maximally 
entangled states, one can find U £ SU(A) such that it 
will never satisfy the preceding constraint. Hence, two 
{SWAP) a gates and single-qubit rotations are not suffi- 
cient to construct an arbitrary U £ SC/(4). 

Spin-based quantum computation by SWAP circuit 
In this section, we will estimate the real time required 
to implement our (SWAP) 01 circuit. We consider both 
the GaAs and Si semiconductor systems, for review see 
Ref.Q. A reasonable effective resonant magnetic field 
is B = lmT, giving a Rabi frequency for an electron in 
GaAs (Si) of approximately 6.2MHz ( 28MHz ), so that 
a single qubit 7r rotation requires approximately 80ns 
(18ns). The effective resonant magnetic field could be 
generated by localized magnetic excitation [3 or through 
g-tensor modulation 0, U- A SWAP gate needs about 
50ps for J w O.lmcV (01, and hence 50ps is the maximum 
time required for performing a (SWAP) a gate. 

From the point of view of operating time (and hence 
qubit storage errors), reducing the number of single qubit 
gates becomes the overwhelming consideration in design- 
ing a circuit. Our circuit is formed by at most three 
(SWAP) 01 gates and six single qubit rotations. Con- 
sidering that single qubit rotations can in principle be 
performed on individual qubits simultaneously, the total 
time to implement a general (SWAP) a circuit is at most 
the time required for three single qubit rotations. The 
optimal CNOT circuit given by Ref.[10( contains three 
CNOT gates and eight single qubit rotations. If we as- 
sume that a CNOT gate takes almost the same time as 
a SWAP gate, we find that the time to implement the 
CNOT circuit is also at most the time to implement three 
single qubit rotations (distinct from the single qubit ro- 
tations of the (SWAP) a circuit). In practice, one would 
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like to implement specific instances of U S SU(4), so 
that gate counts and timings depend on the target op- 
eration U. The self-evident example is that it is best to 
implement U = CNOT with a CNOT circuit, and simi- 
larly U = SWAP with a (SWAP) 01 circuit. The essential 
point is that given the physical means to implement sin- 
gle qubit rotations and (SWAP) a , we can arrive at the 
most efficient decomposition (in terms of operating time 
and gate count) of a target U G SU(4) that might be 
specified using a network of CNOT's, controllcd-phasc 
gates or other logically convenient gates. This optimiza- 
tion is essential for minimizing both storage and gate 
errors introduced by performing the operation U. 

The entanglement power of (SWAP) a gate is strictly 
less than the CNOT gate, but (SWAP)" gate is as ef- 
ficient as CNOT gate (in terms of gate count) in per- 
forming two-qubit operations. From the point of view of 
operating time, we also showed that both (SWAP) a cir- 
cuit and CNOT circuit is the same for a general two-qubit 
unitary operation. So the optimal (SWAP) a circuit is 
as good as the optimal CNOT circuit if the complexity 
of physical implementation of two-qubit gates (CNOT 
gate and (SWAP)" gate) and single-qubit gates are the 
same. However if we directly replace a CNOT gate by 
two (SWAP) 1 / 2 gates and some single qubit gates to 
implement a general two-qubit operation, we need six 
(SWAP) 1 / 2 gates and the operating time will be at least 
doubled compared with the optimal (SWAP)" circuit. 
In this sense the optimizition of (SWAP)" circuit is nec- 
essary, the advantages of the optimization are: first, the 
number of two-qubit gates is reduced by half; second, the 
operating time is reduced to one half to one third due to 
the reduction of the number of single qubit gates. 

Our circuit realizes the three free parameters in 
Eq.lO by adjusting the a in (SWAP)" gate. So 
the SWAP circuit generally involves three different 
(SWAP)" gates with fixed single qubit rotations. This 
is in contrast with the CNOT circuit, in which the two- 
qubit operation, the CNOT gate, is fixed, and the sin- 
gle qubit rotations are tuned according to the free pa- 
rameters in Ea. l|12|l . This is a trade-off in designing a 
circuit, which do we prefer: tuning (SWAP)" gates or 
tuning single qubit rotations? As is generally accepted, 
in solid state with Heisenberg exchange interaction, the 
(SWAP)" gate can be realized by simply controlling the 
interaction time as presented in An arbitrary sin- 
gle qubit rotation is potentially more difficult to realize 
and may take up to 100ns ji^l time to process. For ex- 
ample, one can apply magnetic time varying magnetic 
fields. Since the single qubit rotations are fixed in the 
SWAP circuit of JTBJ, (the single qubit rotations at the 
front and back of the circuit still need be adjusted ac- 
cording to the specified operations), the directions, gra- 
dients, pulse time etc. of the magnetic field can be fixed. 
In this sense, for spin-based quantum computation, tun- 
ing the (SWAP)" gate in our circuit could be advanta- 



geous compared to using a fixed (SWAP)" circuit, such 
as (SWAP) 1 / 2 , and tuning the single qubit rotations. 
Summary and Discussion 

When it comes to solid state implementations, the ex- 
change interaction has emerged as the primary mecha- 
nism for constructing non-local quantum gates, and the 
(SWAP)" gate is the cheapest and the most natural two- 
qubit gate that can be realized using this technology. We 
have shown that simply replacing individual CNOT gates 
with its SWAP circuit is not an efficient implementation 
technique for exchange-interaction based quantum com- 
puting systems. We have presented an alternate opti- 
mization technique and have derived the optimal circuit 
for an arbitary two-qubit unitary operator using SWAP 
gates and single-qubit rotations. 
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Appendix A: Entanglement power of (SWAP)" 
gate As presented in section III, The entanglement power 
of a unitary operator U € C 2 ®C 2 is defined as, 

E p (U) = Average^ i)W2) [£([/|Vi) ® \rh))], (20) 

where the average is over all product states \ipi) ® I ^ 2) € 
C 2 ®C 2 in uniform distribution, see 01 . Zanardi et al\l^ 
showed the following result: 



E P (U) = 1-(C 2 ) 2 h(U), 



fe=0,l 

I k {U) = trT 1+fc ,3 +fc + <^® 2 (T 1+fc> 3 + fc)C^® 2 s T 1 , 3 ) ) 

where the scalar product (A,B) :— tr(AP), T1.3 and 
^2,4 are the transposition operators: T±^ 3 \a,b,c,d) — 
|c, b, a, d), T2 t i\a,b,c,d) = \a, d, c, b). For case U G 
C 2 <g> C 2 , we know that C 2 = 6, and tr7i )3 = trT 2 , 4 = 8. 
With thehelpofr 2 ,4 = (SWAP®SWAP)T l!3 (SWAP® 
SWAP), we know that for U G C 2 ®C 2 , the entanglement 
power of U is 

E p (U) = ^-±[(U® 2 ,T li3 U® 2 T h3 ) 
+ ((SWAP ■ U)® 2 ,T h3 (SWAP ■ U)® 2 T h3 )} (21) 

as presented in section III and in Ref . . 

Our aim is to find the entanglement power of 
(SWAP)". Subsititute U by (SWAP)" in Q, we can 
find the second term in (JSJ) can be calculated as 

(U® 2 ,T h3 U® 2 T h3 ) 
= y + 3 ( e%na + e ~™ a ) + ^{e 2t7ra + e- 2 " a ). (22) 

We can find that SWAP ■ (SWAP)" is similar to 
(SWAP)", we just need to replace e i7r " by -e i7ra in 
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Thus the third term in JSJ is 

((SWAP ■ Uf 2 , T h3 (SWAP ■ U)® 2 T h3 ) 
= II _ 3( e ™« + e - <TO ) + :j( e 2i ™ + e - 2l7rQ ). (23) 

With all of these results we know the entanglement power 
of (SWAP) 01 is 

E p ((SWAP) a ) = 1 - 1 cos(27ra). (24) 

We can find that the entanglement power of (SWAP)" 
is greater than zero except a = and a — 1 which corre- 
sponding to the identity and the full SWAP respectively. 
The \J SWAP has the maximal entanglement power but 
it is less than the CNOT gate. 

Appendix B, Optimal circuit by CNOT gates and sin- 
gle qubit rotations We review here Vidal and Dawson's 
CNOT circuit. We know that the general unitary op- 
erator can be simplified by single-qubit rotations to an 
operator e~ lH which can create arbitrary phases on four 
Bell states. Using CNOT gates, we need to find a circuit 
which can create three independent phases on four Bell 
states; note that an overall phase is not important here. 
It is well known that by CNOT and Hadamard transfor- 
mation, (W ® I)CNOT, we can transform the Bell basis 
to the computational basis |<&+) — * |00), |$~) — > |10), 
|*+) -> |01), |*") -> |11). The Hadamard transforma- 
tion is defined as VT = ^ j ^ ^ /2. The inverse op- 
erator CNOT(kk ® I) transfers the computational basis 
back to Bell basis. Since the operator exp(— iC& z ) ap- 



plies phases |0) 
a circuit, 



|0>, |1) 



s^ll), we can construct 



CNOT(W <g> I)(e- lC2a * <g> e'^ 2a ')CNOT 



-i(l<r z 



)(W ®P)CNOT. (25) 



By this circuit, we can apply arbitrary phases on the four 
Bell states, 



|#+) 
|*+) 
|$-) 
I*") 



»(+Ci+£i+Ca+&0|$+\ 
-i(+Ci-Ci+C3-&)|^+\ 
*(-Ci-fi-C2+6)|$-\ 



-i(-Ci+€i-C 2 -f2) 



(26) 



The parameters Ci?C2;£ij£2 can be chosen so as to re- 
produce exp(— iH) in Eq. \ 1?A thus recovering the result 
given by Vidal and Dawson|lO|. Here we have provided 
a more intuitive understanding of their results. 



The circuit l|25(l is optimal in the sense that the num- 
ber of CNOT gates in it is minimized. Vidal and 
Dawsonflfij. Vatan and Williams ^lj gave two differ- 
ent proofs to show that at least three CNOT gates are 
necessary. Here we provide another proof which com- 
bines both. In Ref.[10(, it is shown that the most general 
circuit involving two CNOT gates can be simplified to 
CNOT(e-^ a * ® e~ l ^)CNOT. In Ref.0, it is shown 
that the simplified circuit is CZ^^"* <g> e^^JCZ, 
where CZ is the controlled-phase gate. Our observation 
here is that these two simplified circuits only have two 
free parameters. With additional single-qubit rotations 
applied before and after the CNOT's, we can have at 
most 12 + 2 = 14 free parameters. But we need at least 
15 since the two-qubit operator acts in a 4-dimensional 
Hilbert space, thus at least three CNOT gates are neces- 
sary. This is a simple proof that the circuit of Eq. [55] is 
optimal. 
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